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Abstract 

The heat transfer processes in a flat plate collector have 
been analysed considering the spectral variation of the radiative 
properties of glass, fluid and collector plate and also the incoming 
solar radiation. The assumption of isothermal ^ass plate has been 
justified by a detailed analysis, A numerical method has been 
developed to calculate the variation of temperature of the glass, 
fluid and collector along the length. The effect of selective 
surface, gDass thickness, glass-collector spacing, fluid (air and 
water) and Reynolds number on the efficiency is discussed. The fact 
thab the flow is not fully developed is shown to be extremely 
important. The advantage in using non-dim ensional numbers is also 
brou^t out. 
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Introduction 


The utilization of unconventional energy sources is going to be 
very important for meeting the future energy demands of our country. 

One of the most important unconventiaial energy source is, of course, 
the solar energy. The energy frcm the sun can be directly converted 
to electricity through solar cells or it can be first converted to 
heat and later to work through heat engines. The most populsir method 
of collectirg solar energy as heat is throu^ the use of flat plate 
collectors. The flat plate collector will continue to dominate the 
field of solar energy collection deices on account of the simplicity 
of its construction and low cost. 

The common flat plate ooilector consists of a blackened metallic 
surface covered with a glass plate (to trap infrared radiation) with 
a transparent fluid flowing betv/een them. Most of the 

work on flat plate collectors was done in the early sixties [see, for 
exaniple, Spanides (l96l)] and interest in this area has increased in 
the last few years. In recent years several modifications have been 
suggested to improve the efficiency of flat plate collectors. They 
include, for example, double ^ass covers, spectrally and directionally 
selective surfaces, radiatively absorbing’ black liquids’ , plastics etc. 
At present not much work has been done to quantify the improvement in 
the efficiency of the collector on accoimt of these modifications. The 
primary reason for this is the insufficient understanding of the basic 
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heat transfer mechanisms. The analysis of the performance of a flat 
plate collector is complicated by the fact tnat the conduction, convection, 
and radiation heat transfer processes interact iti a complex fashion 
(since they occur both in series and in parallel). Major improvements 
in flat plare collector design can be envisaged only after we understand 
the precise roles played by the three heat transfer modes in determining 
the heat loss and gain by the collector. 

The first attempt to analyse the performance of flat plate 
collectors was made by Hottel and Woertz (l942) in which they considered 
the cover plates to be isothermal with the fluid flowing in tubes 
attached belov^ the collector plate. Cobble (l964) determined the 
toaperature field within a transparent solid irradiated by a Planckian, 
or a modified Planckian source. He proposed a simplified model for 
glass which absorbed monochromatic radiation according to Beer's lav/. 
Hov/ever, he considered the other interface to be insulated. Safdari 
(1966) coasidered the transient heat transfer in a transparent material 
in contact with a v/ell stirred fluid placsiin a container. lumsdaine 
(1970) analysed the heat transfer through a semi-transparent medium in 
contact mth a flind which absorbed all incident radiation. He also 
discussed the influence on fluid temperature of the reradiation from 
the solid, the relative importance cf the ambient convective heat 
transfer coefficient and the effect of various physical parameters 
such as solid thickaess, thermal conductivity and the absorption 

coefficient. He did not, however, consider the effect of the variation 
of temperature along the length of the collector. 
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In the present thesis a detailed analysis of heat transfer processes 
in a flat pla’ce collector has been made by taking the following importanx 
aspects into consideration 

1) Isothermal Glass Plate - In most of the previous work on flat 
plate collectors the assumption of the isothermal glass plate has 
been made without a proper justification. In Chapter 2 we shall 
justify tnis assumption by a rigorous analj’-sis based on the work 
of Viskanta and Anderson (l975). 

2) Spectral Uature of Glass •- Ihe spectral nature of glass is 
considered as it will help us in evaluating the merits and dements 
of various types of glass, A simplified model based on the work 
of Cobble (1964) is employed, 

3) Spectral Fature of the Colleotor Plate - It is well known that the 
use of a spectrally selective surface can enhance the efficiency 

of a collector by reducing the infrared radiation emission. An 
analysis taking into account tne spectral variation of the 
collector place emissivity will show the effect of spectrally 
selective surfaces on the performance of flat plate collectors. 

4) Developing Pluid How •- Most of the investigators in this area 
have assumed the flow of the fluid to be fully developed. However, 
for short flat plate collectors the flow will not be fully developed 
and, thus, we have to include the effect of norufully developed region. 

Hon-dimensional Humbers 5?he experimental data that is available 
on flat plate collectors is not easy to evaluate because most of 
them are presented in dimensional foim, The advantages and insist 


5 ) 
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that can be gained by using non-dimensi onal numbers has not been 
fully utilized by workers m this field. As a matter of fact if 
the design of flat plate collectors has to be standardized, as has 
been done with conventional heat exchangers, then there is a need 
to identi:^ the relevant non-dimensional parameters. In order to 
identify the relevant non-dimensional parameters we need to have a 
good understanding of the important heat gam and loss mechanisms. 
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AMIYSIS FOR A RLAO? PLATE COLLECTOR WITH TEMPERATURE 
VARYING ALOIG THE THICKNESS 

We shall restrict our attentioa to a flat plate collector (Pig. 1. l) 
in which a constant property, non-absoihing fluid flows in between the 
collector plateand the glass plate. The collector plate is a diffuse 
emitter of radiation. We assume that the back losses are zero. In 
practice, the back losses can be made to approach a zero value by providing 
a suitable insulation just below the collector plate. The top losses of 
the flat plate collector le dependent on the wind velocity and the conduction 
through the glass, of the energy emitted by the collector plate. Hence, 
the top losses cannot be made to approach a zero value and have to be 
considered in the analysis. The side losses are assumed to be negligible. 

We shall determine the temperature distribution within the glass 
plate. To do so it is neceesaiy to determine the temperatures of the 
fluid and the collector plate. These can be found from the energy 
balances for the fluid and collector plate respectively. We shall, 
therefore, write the equations of the energy balances for the fluid 
and the collector plate before determining the temperature distribution 
within the glass plate. 

An attempt has also been made to determine the variation of the 
fluid, collector plate and glass jlate temperatures along the length 
of the flat plate collector. This is done by dividing the flat plate 
collector into a nuntoer of small elements. The element nearest the 
entrance of the fluid is considered first, Eor the particular element 
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we detenaxiie the temperatui’es of tne glass plate, tne collector plate 
and tne tciiiperature of the fluid ^t che exit of the eleneat. This fluid 
temperature is lov used to determine the temperatures of the glass plare 
aiid tne collector plaie of tho next elemetix. 'le proceed in Idiis manner 
until tae exiL of the fluid is reached. 

In manj’’ collectors, since they are ^lOt very long, tne flow is not 
fully developed and, merefoi'e, one needs to accouit for the variation 
of bne heat transfer coefficient. The variation of the convective heat 
transfer coefficients for a developing laainer or turbulent flow is 
taken from the andysis of Kays (l966). Moreover, the heat transfer 
coefficients at surface 2 of the glass plate and the collector plats 
have been assumed to be the same. 

2.1) Therg/ Balance for the Huid - 

Consider a fluid element Ax at a distance x from the entrance. 
The total convective energy entering the control volume (fig.l.2 ) 


"'^f %f ^^3 (^3 ~ -p 


The total conveccive energy leaving the control volume 


h hf h h to V ''f “ 


Applyiig the conservation of energy we get 

dT 

"f V " ^f ^ ^ ^"3 -"P - ^ 2 ("f " "2^ 


(2. 1) 


I)c[n. (2,1) written in non-i-dimensional fotm is 
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dT, 


dz 


+ 2 


:?u^(x) + -^^(x) 

I 

I 

Fu^(x) T^(x) + nu2(x) T^d) 

Gz 

^f - 2 ' 

Gz 


(2-2) 


where 


2 h^ d 2h d 

f = == VIj, &x ~ Ly/Lf Ft!^ = ~ — , NUg = ,fe=2Re Pr d/L. 


X 


f _ f 

I li 

Pne integrating factor is ezp [ / g{C) and g(c) =2 

0 

The solution to Bin. (2.2) is 

A 

T (x) = T„(o) ezp [-p(x) ] + exp [ ~p(x) 3 [/ P(p) dnJ 


1 u^(?)+Fu 2(?) 

Gz 


(2.3) 


X 


where p(x) = J g(5) d^ 
0 


aiad P(n) = 2 


Fu^(r,) T^(p) + Fu^Cn) 
Gz 



n 

exp 

/ g(OdC 
_ 0 _ 


(2.4) 


The tenperdture of the flmd rit x + Ax is ^iven by 
T„(x + Ax) = T (0) exp [ ~p(x + Ax) J + exp [ -p(z + ax)] 


x+Ax 

J p(p ) dn 

o 


(2.5) 


Subtracbing r-q,n. (2.3) f^^om E^n. (2.5) and expanding p(x + Ax) in a 
Taylor's senes (neglecting second and higher order terms) we get 

T (x + Ax) - T (x) = T (0) [exp { - ^ Ax} - 1 ] exp [ -p(x) ] 


exp [•-p(x)] exp (- ■^ Ax) 


dx 


dx 

X 


x+Ax 


/ p(n) dn + J P(n) dn 

Lo X 


X 

- exp [ -p(x) ] / p(n) dn 


( 2 . 6 ) 
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^ « 1 may replace exp by its series. Neglecting 

dx 

second and hi^er carder terms Eqn. (2.6) reduces to 

T (5 + Ai) = I „(x) - L(0) exp [ -p(x) J ^ Ax 
I i dx 

A 

+ exp [-p(x) ] {1-“T ^ { I ^ 

dx 0 

A 

X 

- exp [ -p(x) J J N(n) dp 


(2.7) 


or 0}^(x + Ax) 


or ^^(x + Ax) 


i) = i (x) [ 1 - ^Ax] + P(x) Ax exp [ ~p(x) ] + O(Ax^) 


dx 


= n(x) 


Nu_(x) + Nu (x) 

[1-2 (-2-ts; ‘icJ 


U . - - (52 


2 C 


m 


■U^(x) ^^(x) + Mig(x) Tg(x) 


Gz 


I 

J Ax 


( 2 . 8 ) 


The above equation can be used to determine (none no ally) the 
fluid tern per aivcce at ar^y point provided the auid, collector and &lass 
temperatures are known at points upstream. Tne coHector and glass 
temperatures can be determined only if tb^ appropriate energy balance 
equations are written for them. This is done next. 

2. ? "^ Energy Balance for the Goliector Plate 

The radiant heat flnx leaving the collector plate, in a 

small frequency interval dv about the frequency v is given by 


® 3 v '^''= "jv 


(2.9) 


IS the radiant heat flrx arriving 


where 


at the collector plate. 
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Si^ialarly, for the radiaat he A flur leaving interface 2 of the 
glass plate, we have 


dv = s e (T ) dv + p H dv + (b) dv 

2v Sv bv 2 gv 2 \) Rv 


( 2 . 10 ) 


wheze is thfe radianc heat flux arriving at interface 2 of the glass 

plate and q (b) is the radiant flux at interface 2 of the glass plete 
Rv 

(see AppSiidix a) . 

¥e assime tnat all the radiant energy leaving a particular element 
of iihe collector plate reaches Lne corresponding element of interface 2 
of the glass plate and vice-versa. iDhen 




dv = , dv and dv = B, dv 


2v 


■2v 


5v 


( 2 . 11 ) 


Erom Eqns. (2.9), (,2. 10) and (2, 1l) we get the following expressions for 


B^ and B„ 
2v 3v 


_ gv bv 2: gv 3^ bv y ^pv 

•“2V "■ 

1—00 

^ ^gv^3v 


( 2 . 12 ) 


®3v ■" 




^ ” ‘^gv ^3v 


(2.13) 


AS the collector plate is insulated, fcne total gain in radiant energy by 
the collector plaLe must be eqiaal to the total energy lost by convection. 


H3 (I3 - V = / -V 


(2. M) 
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^‘5 / 1-V3V i 1-Pgv<^v 


dv 


"^ / 1 - P o_ ^ 

0 "^gv 5v 


The integration over tt'equancy Ccn be resolved into two ranges i.e. , 
the opaque and the transparent. This facilitates the dropping of some of 
the terms, f'or most glasses the transparent region lies between 0-4 and 
4p while the opaque region exists above 4vi. 


or h^(T^-T^) 


/ 

Av 


trans. 


(1 - P ) e e (T ) 

^ Sv* gv bv^ 2 ^ , 

' - “gv “ 3 v 


- / 

Av 




(1 - P ) e-u (l,) 

^ gv' by y 


trans. 


Per ® 

ev 5 V 


/ ^''"^3v^Rv n 

^ / TT ' V - r 

""trans.' ^v 




r '■ "3^ gv''bv''''2^ T r ^3'^^"'’”°gv^^bv^^3^ 

+ J p" -' ^p dv - J -j— 

^'"op. 5v gv 3v 


dv 


t I 

Av 


op. 


'' "■ % ^3v 


dv 


(2.15) 


The last term in Oqn. (2.15) is equal to s^ero as the transmittivoty of 
glass in the opaque regions is zero, is the fractions of the radiant 
energy lying in the transparent regions (o,4u to 4y niost glasses) 
for the ^ass and the collector plates are very small, the first and 
second terms are negligible compared to the other terms. 
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h,(a; -I ) . / 


dv 


^ C e^(5j “• (T,) J 

+ / Sv 3v *- IV 2 ^ by^ 3 

^3v " ^gv '^3v 


(2.16) 


Erom Eqn, (2.16) we do not get the vdli® ol the collector temperature, T^, 

explicitly. In order to do so we make use of tne concept of the radiative 

heat transfer coefficient and replace the last term in Eqn, (2.16) by 

h„ (l- - 0?„) wnere 
3r 2 p" 


^3r " 


e e 
gv 3 V 


£ + S-z *“ £ Err 

gv 3v gv 3v 


op, 


CT (I^ + T^) (T^ + T^) (2.17) 


Eqn. (2. l6) reduces to 
h^ (1^ - I^) = / 


(1 -• ^3y) 


Am ■* - Prz P X 

^"^triTS. ^ 


dv + h^^ (l^ - T^) (2, 18) 


or h (!„ - I,) = / 

P P Av 


trans. 


(l - p ) q (b) 

— ^ ^ dv + h_ (T„-T^)+h_(l_-I ) 

1 - P3^ 3r 2 3^ 3 f 2 


r 


or (aij-ig) = 


hJs„-!P„) + / 


Av, “ ^gv *^3u 

trans. 


"3^ f 2 

2.3) Temperature Distribution \vithia the Glass Elate 


b -‘- 3vbR>) 


(2. 19) 


The general energy equation for steady state and in the 
absence of internal heat generation simplifies to 

|r (9g gp =. 0 


(2.20) 



12 


where ^ and = f dv (2.21) 

In order to reduce Eqn. (2. 2l) into a aon-dnaerxsional form, let 

y = y/b, T = (2.22) 


where denotes the total direct solar neat flux. 


2 - 
d T 

<3y 


qg b 


S 




<iy 


Integrating twice v/ith respect to f we get 


(2.23) 




<is b . 

" K I ^ 

8 “ 


and 


' iV / 

g eo 


(2.24) 


(2.25) 


Eqn. (2,25) specifies the non-dimensional temperature distribution within 
the glass plate, a and 0 are constants which are to be deteimned frem 
the boundary conditions ab the tv/o interfaces of the glass plate. Ihe 
boundary conditions are 


^y ~ ~ ^1 cv^'’“Av ^ 


E 1^, dv 
gv Idv 


op. 


op. 


+ / e [ ©K ~ (l ) ] ^v 

gv bv' V bv 

op. 




op. 


e + E-z "• £„ e,,. 

gv 3v gv 3v 


(2.26) 


dv (2.27) 
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Eqn, (2.27) contains the fluid temperature, the collector plate 
temperature, and the Planck's function of the collector temperature. 
They can. be found from the equations of the energy balance for the 
fluid and the collector plate respectively. 

Using the radiative heat transfer coefficient concept v/e may rev/rite 
Eqns, (2,26) and ( 2 , 27 ) as 


K 


g dy 


y=0 


(h, + h, )(T, 
'‘1 ir' '■ 1 


ni) - I e 1° dv 
^ Av g^ lev 
op. 


- / 


Av 


ir e 1 , dv 
gv Kiv 


op. 


(2.2B) 


y:=D 


(2. 29 ) 


where 


h 


= [ e ] a(l? + T^) (1 1 T ) 

^ gv 1 CO-' 1 

op. 


1 r 


(2, 30) 


Substitution of Eqn, ( 2 . 19 ) in Eqn, (2.29) results in 


K 


g dy 




y=b 


h_ h_ h . 

_ r h + (T - T ) + ^ / 

2 h^^+h^ f 2 ^3 + ^5r ^3v 


(2.31) 


Eqns, (2,28) and ( 2 . 3 l) in non-dime nsional form are 


dlT 


41 = B“i (5^1 - 1 ) - inr , 

^ y =0 8 - iv 


r f ^ / Tt “Sii — ilildy] 

*^3 Av '^d 


op. 


op. 


(2.32) 


and ^ 1 - Bo (T^ „ J ) + ^ J 1 . p 


q_ b h 

^S 3 r 


(l-P,;,,)a-n,,(l) 


^1 


(2.33) 


dv 
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where the modified Biot numbers are defined as 


ii. + h. 




^ =“2 = (‘■a - h^) ^ 

e 5 3r g 


(2.34) 


On substitution of Eqru (2.24) into E^ns. (2.32) and (2.33) 






K T %|* + ^ - Bo (T -1) -■ jr j [ J 

g « y=0 ' ' g ® Av 


~ o 
- 1 

gv loy 


op. 


dv 


+ -2: J ^ 
'^S Au 

op. 


gu Idy 

<ld 


dv ] 


(2.35) 


and 


K I % U 


g 


A A *lq ^ 

+ A = Bo^(T^-Tj + ^ 


( 


3r 


) I 


y=i 


(i-’Pgv)qj^^(i) 


2" f 2" K I 'h-+h_ ^ a! 1 P P 

g„ 3 3r Pg^P 


dv 


^gv"3v 

(2.36) 


where is the total diffuse solar radiation flux. 


Erom Eqn, (2, 25) we may say that 


and 


A A 



q ~ b 

_^S 

K 1 
g “ 


K T 
g ® 


%^y!y=o ° 

% '^y|y=1 + A + 0 


(2.37) 

(2.38) 


On substitution of Eqn. (2,3?) m Eqn, (2.35) and Eqn. (2.38) in 

Eqn. (2.36) we get the following two equations in terms of the constants 


A and C 
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j-j-q I +A=Bo 0- 1} 

g » y=0 g “ 




5nr [ I / irii 


g “ Av, 


op. 


‘S Av 


op. 


'Is 


(Ic- 


9eL + -kV/ 


K T , 

g “ y=1 


V: ^3^3r d. 




1 - P P 
trans. gv 3^ 


Eoam the above two equa cions we get 


A 


[ Bo^Bo^Cl^-l) + ^ ( Bo^(Bo^ / Jj, - qj.^ 


El q, 

J ..E^L-lS^ dv 


E . 1 


Av 


op. 


'Is 


'is 


;,^^dv) -Bo^cg 


h (l-p )q (l) 

, I 3r f 3v -Kv 

+ BOp / Q.T3 ^ylv-i “ h +h 1 1 p p 

tT’ans. 


(Bo^ + Bo^ + Bo^ Bo^) 


and 


(2.39) 


dv (2,40) 


A 

y=l 

dv)}] / 

(2.4l) 
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0 = 1 - 


K 


^ / qg dy|.^ - y - / 

” AV^ 


e X, 

.SI .- JS L a„ 


op. 


q^ - 


'Is Av 


/ ^ - Bo^ - 1) 


op. 


h 


+ 


gr 


(1 - 




dv} 


/ (Bo^+Bo^+Bo^BOg) 


(2.42) 


On substxtutxng for the oonatants from Eqn. (2.4l) and (2.42) respectively 
into Biti. (2. S) we obtain xhe temperature distribution within the glass 
plate. 

Erom the temperat\jre profile within the glass plate shown in 
ilg, (2.1), T'e ftnd that the temperature difference across the glass 
plate is not very large. This is because 1) the conductive resisianoe 
of the glass plate is small compared to the convective resistance at 
the outer surface 2) the contribution of solar absorption to the 
energy balance c£ the glass plate is not substantial and 3) the 
radiative energy exchange between glass and collector are almost equal 
and hence the net radiative transfer from collector to glass is 
negligible. 



GHAHER 3 


ANALYSIS FOR A ILAT HATE CQLLEOTOR WITH TEfflEERATHRE 
CONSTANT ALONG THE THICKNESS 

The results of the previous chapter show that the temperature 
difference across the glass plate is not very large. We may, therefore, 
assume the glass plate to be isothermal. In this chapter we present 
an analysis of the flat plate collector performance assuming an 
isothenneQ. glass plate, 'The other assumptions made are 1iie same as 
that in Chapter 2, One of the major assumptions made in Chapter 2 
was that the fluid flowing between the glass cover and 1he collector 
was transparent to radiation. This assumption is satisfied by air 
but not by water, A one centimeter thick layer of water is opaque to 
all radiation above one micron wavelength while being transparent to 
radiation below one micron wavelength. 

3, 1 Energy Balance for the Glass Hate ; 

The total radiant energy arriving at surface 1 of the glass plate, 
H^, IS the sum of the total direct and diffuse solar radiant energy and 
the radiant energy from the surrounding environment, 

* H, = /° H dv= /” 1 ° dv+ ir f i° dv + /”e (T ) dv (3.1 ) 

. . 1 ' 1v o o o 

The total radiant energy leaving surface ^ of the glass plate, is 

composed of the energy emitted by the glass plate, the fraction of the 
total radiant energy arriving at surface 1 that is reflected, and the 
fraction of the total radiant energy leaving the collector plate that 
is transmitted through the glass plate. 
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B = / B dv = /“ e 6 (T ) dv + / r (yO i° dv 
1 n 1v gv bv ^ ^ gcv^^ ^ lev 


W W ^ 

+ IT J r ,,, 1° dv + / e,^ fT ) r , dv + f t ^ B, dv (3.2) 

gdv 1 dv ^ bv ' gdv ■' gdv 3 v / 


wJiere is the radiant energy leaving the collector plate in a small 
frequency interval dv aroimd the frequency v . The equation for B_,. is 

3v 

V = "bv (V * - SHv 


The first term is the emission from the collector plate and the second 
term represents the fraction of the radiant energy arriving at the 
collector plate that is reflected. is the energy arriving at the 

collector plate. In order to determine we must first determine 
the energy leaving surface 2 of the glass plate. 


®2v = 


e e, (T ) 
gv bv 


+ r, +t (y’) 

gdv 2 V gc V ^ 


0 

1 . + IT 

lev 


t i + t , e (T ) 
gdv Idv gdv bv <» 


(3.4) 


The first term on the r.h.s. of Eqn. (3-4) represents the emission 
from surface 2, the second term represents the energy reflected by 
surface 2, and the rest represent the energy transmitted throu^ the 
glass plabe, denotes the energy from the collector plate that 

arrives at surface 2 of the glass plate. If we make the same assumptions 
regarding the energy arriving at the collector plate and surface 2 of 
the glass plate as in Chapter 2 we have 


2 V 


and H_ = B 
2v 3 V 3 V 


(3.5) 
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Substitution of Eqn. (3.5) into Eqns. (3.3) and (3.4) results in two 
equations in and ^ . Solving these we get the following expressions 
for Bg ^ and B^ ^ 


®2V= [ *gov<'‘'^ "■ ' ^gax. + V'’ "gv°JV 

■*' ^gdv ^3v ~ " ®3vb 

doa Bjy “t®3V % ^^^3^ + (1 - ^goj'^ ^ ^Icy 


■'gdv ldv gdv bv 


Ihe gain in total radiant energy by the ^ass plate must be equal to the 
energy removed by convection. 


hj (Tg - If) + - !__) = (H^ - B,) + (H^ - B^) 


(5.8) 


Where E = j H dv and B = J B dv (5*9) 

^ 0 0 

On substitution of Eqns. (3.l), (3.2), (3.5) > (3*6) and (3*7) into 
Bqn. (3.8) we get [see Bqn. (C.4) Appendix 0 ] 
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1 - Vv 


“ t ^ (1 - e^ ) e [ e, (t ) - e, (T ) ] 

+ J gdV 3V gV^ bv^ oa^ by^ g'' 

° 1 -Vv<’-= 3 v^ 


/ ^ (y')3,?Md'’+ / ire i° dv 

Q gv ^ ^ Icv ^ gv Idv 


"■fc (uO(l~e )c 

+ / .,,SPV— ] g ^. JQ^ dv 


’ - Vv*' -" 3 V> 


“t, (l-e,)e 1 ° 

-Si^ 3\/ gv Idv 

0 1-r, (i~e ) 

gdv 3v' 


(5.10) 


The first term on the right hand side of Eqn. (j.io) represents the 
direct radiant energy exchange between the surrounding environment and 
the glass plate, the second term that between 
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itO'a -Llq_n, ( 3 . 12 ) we 0 ^ deteraiinB tne temperature of the isothermal 
glass plate provided vre knov/ the temperatures of the fluid and the 
collector pic ce, Inese teiii|,eratures are dexerniined from the energy 
balances for the fluid and tne collector plai'^ respectively. The 
energy balance for -che flmd is given by Eqn. (2.8). ¥e, therefore, 
now proceed to 1:116 task of v/mting the energy balance for the collector 
pi ate. 


3 . 2 Energy Balance for the Gollector Hate 


The gain in total radiant energy by the collector plate is, 
in the case of steady state, equal to the heat removed from it by 
convection. 


^ 3 ( 3)3 -Tp = H^-B 3 (3.15) 

Substitution of Eqns. (3.5), (3.6) and (3.7) 1 ito Eqn. (3.13) results 
in {see Eqn. (c.5) Appendiy 0} 


CO £ 


- ’-'f) 


J -SI 




dv 


H- / 

o 


“ ^ r (^ ) - ®v, (^a) 1 

3v gdv ^ b V 0 °'^ by 3 ^ 

*■ ^gdv " ^3v^ 


dv 


e~ t (v') 1. 

+ J .■ 3v ^ cv 

^ ^5v V . v hdy 


I i 5 

CEImIU. 


(3^14^ 

Stiff 

- ♦ U j-tt. 
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The firsTJ term on the rx^t hand sxde of lijn. ( 3 , -]4) represents 
the radxai:ifc exchange of energy'- between the ^ass plaije and the collector 
plate and the second that bet\/een tne environment and the collector 
plate. Tne last two terms represent the solar energy absorbed by rhe 
collector plate, 

Eqn. (3.14-) may also be written ae 


h^CTj - T^ = / 




^ - Vv -"3^} 


dv 


r ^3v ''^gdv f^®bv^^co) “ ®bv^^3^^ 

+ ) ° dv 

^^trans. “ ^gdv " ^3v^ 


^ J ^3v Scv 

trans, ^ ~ ^gdv^^'^ 3 v^ 


dv 


e t X 
3v gdv 1 dv 


+ TT J 

Av 1 - r (1 ~ ^ ) 

trans. gdv ^ 3 v' 


dv 


(3.15) 


The third term and the second term on the right hand side of Eqns. ( 3 . 12) 
and ( 3 . 15 ) respectively are negligible as compared to the other terms. 
Moreover, the solution for the temperatures of the glass plate and the 
collector plate are not easily obtained from the above equations as they 
are non-line ar in nature. We, therefore, linearize both equations 
(using the radiative heat transfer coefficient concept) as a first 
approximation, Eqns. (3.12) and (3.15) now reduce to 
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(h^+hg) 


and 




or !E 

g 


"g - (^ * “2 ’^f) - - V + - V 


+ J ®rr , ( P' ) 1? dp + TT / E i°, dp 

^ gv ' ^ 1cv J gv Idv 


+ 


t_(y') (1 - ^ ) ^.(v’) X 


gcv 


3^)^ 1QV 


Av 


trans. 


1 - Vp ""Sp) 


dv 


•P / ir ''ga !i... l: _L 3 p^ "sp - iigy a„ ( 3 . 1 s) 

Av 1— r ( 1 — e) 

trans, gdv '■ 


^ 5 ^ '*'gov^^ ^ ^ 1 ov 


,) = h„^ (^ -T,) + / 

S3 1-r ( 1 -^) 

trans. gdv^ 5V^ 


dv 


J 

Av 


^3v gdv ^idv 
trans. ^ ~ ^gdv^^ “ 


dv 


( 3 . 17 ) 


[ (h +h, ) T + (h^ + h, ) + / E (u') 1° dv 

\ 1 03 \ 2 3r'^ 3 ^ gv ^ ' lev 


+ ir / 
0 


£ 1 -, dv 

gv Idv 


. / 

Av 


(1 -^^v) ^ 


gov 


3 v/ 


trans. 


1 “ r 


(1 - ^,..) 


dv 


gdv ' ' 3v' 

o 
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Av 


/ , V^»^^'^ 3 p>'^ap 4 w a„] /(h^+b^+h^^^hj^ 


trans 


. ■' " " ^ 3 v^ 


( 3 . 18 ) 
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and 




t (v') 1 


0 

lev 


® 3 = fVf / , , 

^trans. ^ ~ ~ 


d vJ 


<^ + i 

r 3'^ jSldv 1dv I ,, 

J ^ S_^ ^ dvj/(h + h ) 

^\rans. 1 


(5.19) 


Eeam Eqn, (3»19) we determine the collector temperature, T^, in terms 
of the glass plate temperature, T . Substitution of this value of the 

O 

collector temperature, f , m terms of the glass plate temperature, I , 

in Eqn. (3. 18) results in an expression in the glass temperature and 

tne fluid temperature, I , Using the value of the glass temperature, T , 

1 g' 

obtained from this expression we determine the collector temperature 
from Ikpa. (3.15), This collector temperature is now used in Bq.a. (3,12) 
to determine a newr glass temperature. We now go bach to Eqn. (3.15) and 
determine a nevf collector temperature. iThe process is repeated till a 
reasonable convergence of the glass plate and collector plate temperatures 
are obtained. 


Once the temperature of the fluid at tne exit is known, ae 
can determine the efficiency of the flat plate collector from the 
equation given below 


(T 


- T ) 
f,oat f,in^ 


rj 

_ O 


dv + n f 1 °^^ d'^ X 1 


Efficiency = 


0 

1^ 

lev 


(3.20) 
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or on making the appropriate transformations 


Efficiency 


1 (-1) X (. 

2 


Gz 


“) X (T 


/ 1*^ dV + TT / l'!^-... "iV 

■' 1cv > 


0 

'Idv 


f ,out 


f- ) 

f ,ia 


(3.2l) 


3.3 Radiativsly Absorbing Fluid 

The analysis given previously is not valid for an absorbing 
fluid like water. The absorption coefficient in the region less •cnsn 1y 
is negligible as shown in lig. (3.2 ) which is taken from Goody (l964). 

It is a reasonable assumption that water transmits all radiation at 
wavelengths below 1 y and it absorbs all radiation at wavelengths 
greater than 1 y . Hence, as far as tne emission of radiation is 
concerned, water is assumed to behave like a blackbody. We also 
assume that the water is isothermal (at T^) for the sake of simplicity 
of analysis, we now proceed to the task of v/riting the equations ox 
the energy balances for the glass plate, the collector plate, and 
the fluid respectively based on this simplified model. 

a) Energy balance for the glass pLate ; 

The energy absorbed by the glass plate must be equal to the 
sum of the energy emitted by it and the convective energsr removed 

from the two surfaces. 



(a I 
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4 X 

aTp 


2 e 


g 








{5. 22) 
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b) Energy balance for the collector plate 

Ihe collector plate absoros rhe solar energy transmitted first 
through the glass plate and rhen through water. It also absorbs energy 
emitted by the water. The loss is through emission and by convection 
at its surface. An energy balance, therefore, results in the following 
expression 


I e 1 dv + 

0.4p 


it/* e 1° dv+e oT^ 
0.4P Idv 3 f 


3 


oT: + h,(T,-Tj 


3 '‘■'3 f 

(3.23) 


c) Energy balance for the fluid 

An energy balance for a fluid element of length dx results in 




4y - 0 

+ { / t 1 ° dv + ir / t 1 dV } dx 

IP 1U 


d T^ 

■f 'f ''ptf '~f ' IT 


2 otI dx + d 0,^ (T^ + -~r~) dx 




(3.24) 


dT „ 4 

h " V sr = ■'a ‘"g - - "p * =g '“'g 
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+ e o(T^ - 'F) + / t 1 hv 


3 

4V 


+ / TT t , 1^ , dv 

gdv idv 
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(3.25) 
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Eqn, (3* 25) is non-dimensional foim is 


2Eu2(x) 


dx 


2Fu^(x) ^ 
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2 f m^ m^ 'k . ^ 
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{ ~ } + Tr_ T 02 / 
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(3.26) 


V/e integrate Eqn. (3.26) to get 5^(x). Using the method given 
in Section (2.1) we get T^(x + Ax) as 


!I}^(x + Ax) = !r^(x) + { i^(x) + ®^(x) - 2l^(x) } ax 


2011 (x) 


R„ ^ (ll.(x) - Tt(x)} + Ej {(T3 (x) - ^^(x)} 


Gz I g ' g 


^<18 ^ . ^ov . . ^^d .. Ndv, 

+ ~ ^ n TT T (l!7, / ^ 


Gz “gov Ig Uz gdv 1 , 


■d 


T Gz 

m - >4 

(3.27) 

In Eqns. (3.22), (3.23) and (3.27) we have three non-linear 
expisssions in the ^ass plate, collector plate and the fluid temperatures. 
Ihese three unknown temperatures may be found using the procedure 


explained earlier. 



GHAPCBR 4 


RESULTS km DISCUSSIONS 


The equations presented in the last chapter "-ere solved 

numerically using the following input data s 

1 ) The spectral variation of aosorption oceffioient and reflectivity 
for different types of glass were taken from Gleek (1959)» 

2 ) The heat transfer coefficient on the outer surface of glass was 
taken to he equal to ( 5.7 + 5»8 V , ) W/m^ °C as suggested hy 
Duffle and Beckman (l974). 

3 ) The variation of Nusselt number on the inner surface of glass 
and the collector was taken from the solution for developing 
flov; given hy Kays (1966). The expression for Nusselt number 
IS given as 


= 


1 exp (- 

2 I eicF J/Ez) 


where \ and & 
n n 


are eigenvalues and eigenfunctions respectively. These are 
tabulated in Kays (1966). In the case of laminar flow, where 
the fully developed Nusselt number is not a function of Reynolds 
number, the above expression can be used directly. For turbulent 
flor/, however, we have to include the effect of Reynolds number 
on the fully developed Nusselt number. This was incorporated by 
multiplying the above expression by (^*)*^*^ where Re-*^' is the 
reference Reynolds number for which the eigenvalues and functions 
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•were available. 

4) The spectral variation of incoming direct solar radiation were taken 
from Coulson (1975) and that for diffuse radiation from Probert and 
Hub (1968). 

Prom Pig. (4*1) and Pig. (4 •2) we see that the efficiency of the 
collector increases continuously vnth the Reynolds number. This is 
because the outlet ten^erature (and hence the glass and collector 
temperatures) decrease as the Reynolds number increases and, hence, 
the convective and radiative losses to the ambient also decrease. At 

I 

veiy large Reynolds numbers the efficiency asymptotically approaches 
a value around 75?^ f oi: a water heater and around 80/^ for an air heater. 

In this limit the convective and radiative losses are small and the 
main losses are the optical losses due to reflection and absorption 
in glass. 

Pigs. (4.5) and (4.4) show the variation of the collector, 
glass and fluid temperatures along the le ngtn of the collector for 
air ai^ water respectively. We see that in the v/ater heater, the 
collector temperature is greater than the fluid temperature which is 
greater than the glass temperature. In air heaters hoth the collected 
and glass ten^eratures are well above the fluid temperature. The high 
^ass temperature implies higher heat losses to the ambient. Hence, 
air heaters are much less efficient than water heaters. In water 
heaters the temperature difference between the glass, collector and 
fluid IS very small (around 2°C) while in air heaters the ten^erature 
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differences are of the order of 50° C . This is primarily due to the 
lOT/er heat transfer coefficient of air as ooiiipared to water. 

In Figs. ( 4 . 5 ) ard ( 4 . 6 ) the increase of the temperature of the 

fluid (i.e. , T^ - T^ la plotted against the Reynolds number. 

We see that as the Reynolds number increases the outlet temperature 

decreases rapidly. As a matter of fact for a given application 

(T - T„ ) will be speciiied and, hence, the Reynolds number 

^ f,out f,in 

is fixed. Further discussions on this point will be taken up later. 

In Pigs. ( 4 * 7 ) and (4-8) the effect of emissivity of the collector 
on the efficiency for both air and water heaters are shown. The 
efficiency decreases rapidly with the decrease in emissivity but 
interestingly the efficiency does not become zero as the emissivity 
of the collector goes to zero. This is because even when the collector 
reflects all the radiation, the water is able to absorb some of the 
solar radiation directly and in the case of an aix heater some heat 
IS transferxedfrom glass (v7hioh absorbs solar radiation) to air. 

In Pig. ( 4 * 9 ) the effect of the angle of inoidenoe of direct solar 
radiabion on the efficiency is shown. We see that when the angle is 
less than 45® (y'^ 0*7) the efficiency does not decrease substantially. 
At angles greater than 45° the decrease in efficiency is very large and 
IS on account of the rapid decrease in the transmittivity . A plot of 
the transmittivity (for direct radiation) with ii'(= cos 0*) is shown 
in Pig. (4 .is) and it is clearly seen that the trend shown here is 
the same as in Pig, (4.9). 
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In Fig, (4.10) the effect of acfoient p/ind velocity on the 
efficiency of an air heatt,r is sho'vn. We seo that as the amhient 
wind velocity increases the efficiency decreases rapidly. This is 
■because the most dominant resistance to heat loss from the fluid is 
the convective resistance at the outer surface of the glass plate. 

As this resistance decreeses (with increase in vVind velocity) the 
heat losses increase rapidly. 

In Fig. (4.11) ■ne see the effect of glass thickness on efficiency. 
As the glass thickness increases the efficiency decreases rapidly on 
account of the decrease in transmittance of the solar radiation. One 
may "be tempted to think that as the glass thickness increases the heat 
losses to the am'bient should decrease on account of a greater conductive 
resistance. But as pointed out at the conclusion of Chapter 2 the 
conductive resistance of glass is negligi'ble as compared to the 
convective resistance at the outer face of glass and, hence, an 
increase in glass thickness does not contri'bute to a decrease in heat 
losses "but only to a decrease in transmittance of the solar radiation. 

In Fig. (4.12) the spectral variation of absorption coefficient 
and reflectivity is shown. Glass A absorbs more of the solar rahiation 
(in the region O.4 to 2 .7 p) than Glass B. At higher Reynolds number 
the energy absorbed by the glass is transferred to the fluid more 
readily than at Imr Reynolds number. Hence at high Reynolds numoer 
the Glass A has a higher efficiency than Glass B. 
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In Figs. ( 4 * 13 ) and ( 4 . 14 ) we show the effect of the depth of 

flow on efficiency. As the ratio d/h increases, the efficiency decreases 

on account of decreasing heat transfer coefficient. In the case of air 

the flow IS never fully-developed. In the case of water, the flow is 

(3. ci 

fully developed of O.Ol . As ^ is increased from 0.01 to 0.02, the 

heat transfer coefficient decreases and, hence, the efficiency decreases. 
Increas ing Yj' from 0.C2 to 0,03 does not result in a substantial change in 
effieionoy because of t?<o counteracting effects. As the ratio is 
increased the fully developed heat transfer coefficient decreases but 
at the same time the region of developing flow is more in the case of 
~ = 0,03 than ~ “ 0,02. It is well known that in the regions of 
developing flav the heat transfer coefficients are hi^er than the 
developed region. Hence, the two effects cancel each other 

resulting in no substantial increase in efficiency. 

Figs. (4.15)> ( 4 * 16 ) and (4*17) show some characteristics of 
selective surfaces . The use of selective surfaces is usually resorted 
to in order to increase the efficiency of flat pleto collectors. In 
Pig, ( 4 . 15 ) we conpare the efficiency of the air heater for selective 
surfaces with various cut-off frequencies. We see thab as the cut-off 
frequency is increased from 1 y the efficiency increases. The increase 
in efficiency becomes smaller as we approach 4 v • This is because most 
of the solar radiation lies in the region 0.4y - 4y and hence a selective 
surface which has a high absorptivity from 0 to 4v low absorptivity 

above 4 m (la which region the emission from the collector lies) is the 
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best choice. It is interesting' to note that a collector with emissivitjr 
eq.tial to one can have an efficiency greater or less than a collector with 
a selective surface. At very high Reynolds numbers the collector 
temperature is low and hence the radiative losses are also small. In 
this case a collector with emissivity 6Q.ual to one is better than a 
selective surface because it absorbs all the solar radiation while the 
selective surface absorbs only 90^ of it. In the low Reynolds number 
region, the collector temperature is large and, hence, the heat losses 
are substantial. Therefore, the selective surface (which has a low 
emissivity in the region greater than 4p) ^.s superior to a black surface 
under these conditions. Prom the above discussion we see that at low 
Reynolds number (i.e., high outlet temperatures) we need a selective 
surface with a out-off wavelength of 4v • selective surfaces 

developed so far have a cut-off wavelength around 1 y (e.g. Silicon Oxide 
on Aluminium has a cut-off wavelength of about 1 .6 y) . The above analysis 
thus shows that selective surfaces with higher cut-off wavelengths must 
be developed. 

In Pig.(4.l6)we show a comparison between a selective surface and 
a gray diffuse surface having an emissivity of 0.9. The gray diffuse 
surface can have an efficiency lower or higher than a selective surface 
having a cut-off wavelength of 1 .Oy or 1 .5y depending on whether the 
Reynolds number is low or high. The reason for this behaviour has heen 
explained above. However, a selective surface with a aut-off wavelength 
of 4y has a higher efficiency than that of a gray diffuse surface having 
an emissivity equal to 0.9. This is because, even though the solar radiation 
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absorbed by the two surfaces are equal, the radiative heat losses 
using a gray diffuse surface of emissivity 0,9 is greater than that when 
using a selective surface having a cut-off wavelergth equal to 4 p • 

In Fig. ( 4 . 17 ) a comparison between a black collector and a selective 
surface with a cut-off wavelength of 1 ^ is shown. ¥e see that the 
black surface has a greater efficiency than a selective surface. Increasing 
the cut-off wavelength beyond 1 p will not bring about any change in the 
results as water completely absorbs all the radiant energy lying 
above 1y, 

In Fig. ( 4 . 19 ) we compare the present method with Hottel's method 
as presented in fuffie and Beckman (l974) . We see that Hottel's 

method predicts a lower efficiency than the present method fcr all 
Reynolds numbers, This is primarily due to the fact that the flow 
is not fully developed and hence the local heat transfer coefficient 
is a strong function of the distance along the collector. In 
Hottel's method one is forced to choose a mean heat transfer coefficient 
and this is much lower than the local values near the entrance and, 
hence, the efficiency predicted is lov/er. As the i^ynolds number 
increases the heat transfer coefficient also increases and hence the 


discrepancy is larger, 
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appmdix a 


Por the flat plate collector the emission of glass is 
negligible compared to the solar energy absorption by glass in 
the Visible and near-infra red regions. We can, therefore, 
assume it to be a cold medium. The radiative flux for a t rani us cent 
cold medium la given by Viskanta and Anderson (l975) 


^^o^ ^ ^ ^ (-'C/wo')] 1^0 ""lev V 


+ 2 [ 1°^^ / exp V' dv'/3y(<°,p) 
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K,. + K 


+ igdv ^ dp’/ev(<°,p) ] 


K° - K 




0 

1 2k K 

-2[i°^^ / X^^(p')p2v(p)®^p(-“i; du'/e^{K°,y) 


IC° „ K 


^2dv ^ ^ dp’/e^ (%,p) 3 (a.1) 
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vshere 


= 1 - P2v TT^ 


= 1 - Av)^ (i - v'^) 


(a. 2) 
(A.3) 
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0 o 

and y^ IS the value of y correspoading to y* = y^ xn (a.3)< 

!Ih.e first term on the right hand side of Ejn. (a. l) denotes 
the heat flux due to the collimated beam incident on interface 1 
after talcing into account multiple internal reflection and attenuation, 
!Ilie first term within the parentheses of the second term represents 
the heat flux due to the diffuse ccmponent incident on interface 1 
and the other term takes into account the reflection at interface 1 
of the diffuse component of radiation incident at interface 2. Ihe 
last two terms can be explained similarly. These are, however, negative 
quantities as they are backward and directed. 

Per the case of a flat plate collector, only the first term 
and the first term within the parentheses of the second term on the 
n^t hand side of Sin, (A, i) are of importance as all the other 
terms are negligible as compared to these two. The intensity of 
radiation at interface 2 of the ^ass plate is diffuse in nature 
and does not contain a collimated components Moreover, the radiant 
energy from the collector lying in the regions in which glass is 
transparent to radiation is negligible. Hence, the contribution 
from the collector plate to the radiative heat flux within the glass 
plate Will be negligible as compared to that frem the solar energy, 

Eqh, (A,i), therefore, reduces to 



58 
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lev 


exp(- -^)/6 (<, 
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+ 21 
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Idv 


/ 


T^^(v') exp (-Kyy) y' dy'/g^ (k°, y) 


(A.5) 


Ejn, (a. 5) gives the radiative heat flux within the glass plate in a 
small frequency interval around the frequency v • order to obtain 
the total radiative heat flux within the glass plate, one has to 
integrate Bjeu (A.5) over the entire range of frequency. 
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APPENDIX B 


The reflectivity in the transparent region p see Siegel and 
Howell (1972)J niay he expressed as 


(uf) - Q) 

^ 2 2 / V 

^ sin (0» + e) 


■j + OQS^ (Q* + Q) 
cos^ (0' - 0) 


We also have the following relations s 

Py = Sin 0*/sin 61 y’ = oos 0', and y = cos 9 


(b*1) 


(B.2) 


Prom these relations we get 


tan 0' = {1 - y*^}'Vv’ and tan 6 “ {1 - y’^>*V{’4 - (1 - y'^)}*^ 

(B.4) 

Also 


or 


and 


or 


oos (o' + 0) _ 1 <“ tan 6* tan 8 
cos (e’ - 0) "" 1 + tan 0* tan 8 


oos (9' + 9) 
cos (0' - 0) 


P {n^ - (1 - y«^)}*^ - (1 

y» {n^ - (1 - y*^)}*^ + (1 



sin (0^-0) ^ 1 - cot 9« tan 9 
sin (e* + 0} "" 1 + cot 0* tan 0 


sin (9' ~ o) 
sm (0’ + 0) 


{n^^- (1 - y«^3}'^ -y» 

{n^ - (1 - + V 


(B.5) 


(b *6) 


Substitution of Eqns. (B.5)> (B«5) and (B.6) mto Eqy. (b.i) results in 
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P^(y') = 1 - 2y' {nj - (1 ~ 




y' + {r^i - (1 - y^^)}’^ 


)' 


+ (■ 


n , 


y' nj + {4-(l 


-u'")} 


2^..5' 


(B.7) 


and p (y) = 1 - 2n, y {1 - n^, (1 - y^)}*^ 


( 

L % V. 


+ {1 - (1 - y )} 




+ (- 


y + {1 


_J )2 

- n (1 - y )} ■" J 


(E.8) 


The trans mitt 1 VI ties are then given hy 


T^(y') = 1 - p^(y0 (B.9) 

and T^(y) = 1 - P^Cy) (B.10) 

The effective directional reflectance r^(y') and transmittance, ty(y' ) 
are given helow [[^Viskanta ( 1975)1] 


3 ^^(y') = ■*" exp(-2k^h/y)/B^(k^ h,y) (b.11) 


and ty(y') = '^•^^(w') (b.12) 

The hemispherical reflectance and transmittance for diffuse radiation 
incident on the semitransparent sheet is given by 

1 

r = 2 / r (y') y’ dy' (B»13) 

V i V 

1 

ty = 2 / ty(y') y’ dy» 


(B.14) 
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APPENDIX C 


l) Energy balance for the glass plate . 


The expressions for and B^^ are from Chapter 3 


B.,^, ss e e, (T ) + r (p') i^ + it r i^ + e (T ) r 
I'' gv gov'-^ ^ ^10V gdv idv bv^ g^ gdv 


gdv 

^gdv ®3v 


^1V “ ^lov ^idv ■*■ 


(3-2) 

(3.1) 


®2V“ ■*^^gcv^ ^ ^ ^Icv'*'^ ^gdv^ldv’’’ \dv ^gv %v ^^g^ 

Vv "3V 


®3v = {"3v %v^^3^ ^ *■ "3v^ "bv^V ^ ^Icv 

^ ^ Vv ^Idv ^ Vv %v 1~W-^3v^^ 


Also 


H„.. = B,.. and = B, 


2v " 3v 3v 2v 

On sijbstitubion of Eqn. (3.5) into Bjn. (3.2) we have 

^gcv^'^'^ ^1cv ^ ^g4v ^Ihv ^gdv%v ^^g^ ^gv %v ^^g^ 


1 - 1'^'’ 




gdv^ ^3v 

+ e e. ( T ") 

. gdv 3v ^vy s: 

^ “ ^gdv^^~^3v^ 


(3-5) 


(g.i) 
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ilso 
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w to 
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+e,_ (Tj ] X 
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^ 1 . r {1 . e ) ^ I 


r 0 r ^■' “ ~ ^ 

+ ^ ^ ] dv 


gdV l(3y' , . 

' gdv^' 3V'' 


^ "^gv f^®bv^^»^“®bv^V 

l'^ ^ 




dv 


° ■'-Vv<l-Sv> 


+ j” V'-^'' - ^gy 
Ejn. (c,3) inay be further reduced to 


e, (T )] 
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(C.3) 


05 CO 
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' - Vv - ^3~) 


d V 


or (H^-B,).(H^^^) = / ] do 
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"■ ^gd V ^ ^ 3v^ 


2) Eaergy Balance fot~ the collector plate 
Item Eqas, (3.5) and ( 3 . 7 ) we have 
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or 
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